Analisi dimensionale: due interessanti 

applicazioni 

Germano D'Abramo 

Istituto di Astrofisica Spaziale e Fisica Cosmica, 
Area di Ricerca CNR Tor Vergata, Roma, Italy 
E-mail: dabramoOrm . iasf . cnr . it 

5 Giugno 2003 



1 Introduzione 

Malgrado non sia molto conosciuta persino fra gli addetti ai lavori, V analisi 
dimensionale costituisce uno degli strumenti piu potenti e intriganti delle 
scienze fisiche. 

Ma in cosa consiste? Brevemente, dato un sistema fisico complesso, come 
ad esempio la formazione di un'onda d'urto nucleare o la creazione di un 
cratere d'impatto, e definite tutte le variabili e le costanti fisiche che sembrano 
descrivere propriamente e completamente il processo, spesso risulta possibile 
scrivere in maniera univoca delle leggi matematiche (solitamente leggi di 
potenza) che forniscono la giusta dipendenza funzionale (dal punto di vista 
fisico) di una di queste variabili rispetto a tutte le altre. E questo puo 
essere fatto semplicemente prendendo in considerazione le dimensioni fisiche 
fondamentali delle variabili e costanti fisiche coinvolte: lunghezza, massa, 
tempo, carica elettrica (in alcuni casi specifici anche la temperatura) e le 
loro combinazioni. 

Un'applicazione immediata. Supponiamo di non sapere molto di cine- 
matica e di voler ottenere la relazione matematica che lega l'accelerazione 
centripeta, a, la frequenza di rotazione, u 1 , e la distanza dal centro, r, di 
un corpo che ruota a velocita uniforme lungo una traiettoria perfettamente 



^Xa frequenza di rotazione e defmita come 2ir /T ', dove T e il periodo di rotazione. 
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circolare. Euristicamente parlando (e anche in base alia nostra esperienza 
diretta), e ragionevole immaginare che a dipenda proprio solo da u e r. 
Consideriamo allora le dimensioni fisiche di a,u e r: 

[a] = L-T~\ 

M = T~\ (1) 
[r] = L. 

II segno [x] sta per "le dimensioni fisiche di x" e L sta per unita di 
lunghezza e T per unita di tempo. 

Ora, l'unico modo di ottenere le dimensioni fisiche di a e di combinare uj 
e r nella maniera che segue: 



Va notato che la quantita 



u 2 r. (2) 



(3) 



b adimensionale e che quindi risulta invariante per cambiamento di unita di 
misura, ad esempio passando da metri a chilometri, da secondi a ore. 
Quindi, abbiamo che la e una costante 

k, (4) 



u 2 r 



cioe 



a = kto 2 r, (5) 

che e proprio la dipendenza funzionale che si otterrebbe sviluppando matem- 
aticamente la teoria cinematica del moto circolare uniforme. 

Ovviamente, l'analisi dimensionale puo fornire la dipendenza funzionale 
tra le variabili ma non e capace di definire le costanti moltiplicative adimen- 
sionali, la costante k nel nostro caso. Per queste sono indispensabili o una 
vera teoria fisica o un confronto diretto fra le relazioni ottenute e i dati speri- 
mentali. Va detto inoltre, come vedremo in seguito, che tra i limiti di questo 
metodo di studio c'e che non sempre e possibile definire univocamente la re- 
lazione funzionale tra le variabili coinvolte e anche in questo caso h necessario, 
in mancanza di una vera teoria, un confronto diretto con i dati sperimentali. 
Per inciso, nel semplice caso analizzato sopra vale esattamente a = u 2 r (cioe 
la costante adimensionale k e proprio pari a 1). 
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Piu in generale, supponiamo di avere una relazione funzionale tra una 
quantita a, che deve essere determinata in un esperimento, e un insieme di 
parametri (ai, a n _i), che sono sotto il controllo sperimentale. Sia a che 
l'insieme (ai, a n _i) possiedono dimensioni fisiche. La forma generale della 
relazione puo essere la seguente 



Ora, se il minimo numero di parametri nell'msieme (a,a\, ...,a n _i) con i 
quali posso ottenere tutte le dimensioni fisiche presenti nel mio problema e 
k < n, il Teorema II di Buckingham pQ, che rappresenta il risultato centrale 
dell'analisi dimensionale, afferma che il sistema fisico alio studio dipende da 
un insieme di n — k parametri adimensionali (III, n n _/-). Questi parametri 
sono combinazioni adimensionali indipendenti delle quantita (a, ai, ...,a„_i); 
essi sono indipendenti nel senso che, per ogni i, non e possibile ottenere LTj 
da una combinazione di altri lis. 

In questo modo, l'eq. © puo essere scritta come 



dove <3> e una funzione adimensionale generale. Ancora, per ottenere la forma 
esatta di $ abbiamo bisogno di una vera teoria del fenomeno fisico alio studio, 
o un confronto diretto con i dati sperimantali. 

In questo articolo applichiamo le procedure dell'analisi dimensionale a 
due casi molto diversi fra loro (e con finalita diverse), ma entrambi piuttosto 
interessanti. 

Nel primo caso (sezione 2) si usa l'analisi dimensionale non per derivare 
una relazione funzionale tra quantita fisiche, ma per capire come cambia 
la scala temporale dei fenomeni fisici al cambiare del valore numerico dei 
parametri importanti, quali la gravita superficiale di un pianeta. 

Nel secondo caso (sezione 3), invece, utilizziamo l'analisi dimensionale 
nella maniera piu canonica per ottenere una relazione matematica che leghi 
l'ampiezza di un onda marina alia sua velocita. Applicazione immediata: la 
fisica degli tsunami 2 . 

2 Onde anomale particolarmente distruttive. Letteralmente onde di porto, dal giappo- 
nense tsu (porto) e nami (onda). 




(6) 




(7) 
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2 Moviola sulla Luna 



Avete mai visto, proiettati a velocita doppia, i filmati degli astronauti die si 
muovono sulla superficie lunare? Se provate a farlo vi accorgerete sicuramente 
di una caratteristica stupefacente: sembrano girati sulla Terra. I movimenti 
degli astronauti, i sobbalzi del rover motorizzato etc. sembrano proprio quelli 
naturali che si avrebbero se le stesse operazioni venissero compiute sulla 
Terra. 

E questo un ulteriore elemento a favore della famigerata teoria secondo 
la quale la NASA non avrebbe mai portato l'uomo sulla Luna? 

Molto probabilmente no. In ogni caso qui di seguito mostriamo come 
questo curioso comportamento sia legato semplicemente alia gravita superfi- 
ciale dei pianeti. 

In generale, ogni fenomeno dinamico la cui evoluzione nel tempo dipenda 
anche dalla gravita di un pianeta (ad esempio, un essere umano che cammina 
e salta, un palazzo che crolla sotto il proprio peso, una slavina o una colata di 
magma che scivolano a valle) ha un tempo caratteristico che dipende proprio 
dalla gravita superficiale. II tempo caratteristico non e una costante fisica o 
un numero ben preciso, ma e un concetto generale. Per dare un idea, il tempo 
caratteristico della vita umana e circa 80 anni, mentre quello della dinamica 
orbitale dei pianeti interni del Sistema Solare e dell'ordine dell'anno (periodo 
di rotazione della Terra intorno al Sole). 

Ad esempio, una delle piu semplici equazioni della cinematica ci dice che 
una biglia che cada da ferma in un campo gravitazionale costate g percorre 
nel tempo t uno spazio pari a 



Invertendo l'equazione (jHJ), si ottiene il tempo in funzione dello spazio 
percorso 



L'equazione © ci mostra proprio la dipendenza del tempo caratteristico 
dalla costante gravitazionale g. Infatti, piu grande e g e minore e il tempo 
impiegato dalla biglia per percorrere in caduta libera lo spazio x. 

In generale, quindi, i tempi caratteristici con cui si verificano i fenomeni 
dinamici sulle superfici planetarie con gravita superficiali g\ e g<± scalano come 



x = 




(8) 




(9) 
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h [g^X ■ v rp rp 

T 2 = fe- cloe T ' = 

Nel caso specifico dei filmati lunari il tempo caratteristico lunare T 5 e pari 

a 

T, = T */f = r «/il- 2 ' 45r «' < n > 

dove T®, (7® e so no rispettivamente il tempo caratteristico terrestre, la 
costante di gravita superficiale terrestre e quella lunare (ricordiamo che la 
gravita superficiale lunare e circa 1/6 di quella terrestre). 

Cioe, Tj e poco piu di due volte maggiore del tempo caratteristico dei 
fenomeni dinamici sulla superficie terrestre. E quindi, se aumentiamo di due 
volte la velocita di proiezione dei filmati lunari otteniamo proprio l'impressione 
che i fenomeni dinamici si svolgano con la stessa tempistica dei fenomeni ter- 
restri! Ricordiamo pero che la relazione (fTTj) vale solo per i fenomeni che sono 
soggetti alia forza di gravita: se fosse possibile intravvederle, noteremmo che 
le palpebre degli astronauti battono ad una velocita innaturalmente doppia! 

3 Cavalcare lo tsunami 

In questa sezione, sempre grazie all'ausilio dell'analisi dimensionale, deriver- 
emo la relazione che lega l'ampiezza di un onda marina alia sua velocita. 

Attraverso questa relazione si puo comprendere per linee generali il motivo 
fisico della formazione degli tsunami, ovvero delle tanto temute onde anomale 
che possono raggiungere altezze ragguardevoli in prossimita della riva ed 
essere altamente distruttive. 

Cerchiamo in primo luogo di ottenere l'energia totale che puo essere tras- 
portata da un'onda marina. Come sara chiaro nello sviluppo successivo, 
l'uso dell'analisi dimensionale ci esime dalla necessita di rappresentare (e 
semplificare) Fonda marina come sinusoidale piana, cioe periodica, il cui 
profilo segue perfettamente la funzione trigonometrica seno (o coseno) e il 
cui fronte e una lina retta infinita. 

Visto che abbiamo a che fare con un fenomeno fisico che coinvolge l'acqua 
del mare, cominciamo ad elencare tutte le quantita fisiche che euristicamente 
ci aspettiamo possano entrare in gioco nella descrizione del fenomeno alio 
studio, insieme con le loro unit a di misura. 
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Figure 1: Rappresentazione schematica di un'onda, con le grandezze fisiche 
die la definsicono (per ovvi motivi non e rappresentata la densita p). 

Immaginiamo di essere interessati per prima cosa aU'energia cinetica tras- 
portata dall'onda per unita di superficie attraversata (che chiamiamo E c /S). 
Se parliamo di energia cinetica non possiamo allora non considerare la ve- 
locity dell'onda, v, che sicuramente influisce, cosi come la sua ampiezza, A, e 
la densita del mezzo in cui l'onda si propaga, p, che fanno le veci della massa 
nella classica relazione dell'energia cinetica E = |mw 2 . E ovvio che in questo 
caso non possiamo utilizzare m perche non stiamo parlando di un oggetto 
dai contorni fisici ben defmiti. 

y, P, v, A, A. (12) 

La lista (|T2Jl raccoglie quindi tutte le quantita fisiche considerate impor- 
tanti. In realta compare anche A, che e la lunghezza d'onda dell'onda marina. 
Essa e una caratteristica geometrica dell'onda che non abbiamo ragione di 
escludere. 

Per la derivazione di E c /S non abbiamo preso in considerazione la costante 
di accelerazione gravitazionale g. Ci aspettiamo infatti che per E c /S la 
gravita non sia rilevante e che l'energia cinetica sia determinata una volta 
fissate A, p, v e A. Prenderemo in considerazione g piu avanti quando de- 
riveremo l'energia potenziale (gravitazionale) dell'onda. 

Veniamo ora alle unita di misura (dimensioni) dei parametri in (|12j) : 
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[E c /S] = M-T-\ 
[p] = M-L~\ 

[v] = L-T-\ (13) 
[A] = L, 
[A] = L. 

dove M sta per unita di massa, T unita di tempo e L unita di lunghezza. 

Dall'elenco (fT3|) si puo notare come con i parametri A, p e v e possibile 
costruire, attraverso le giuste combinazioni, le dimensioni di tutti gli altri. 
Essi sono cioe i parametri dimensionali indipendenti del problema secondo la 
teoria dell'analisi dimensionale di Buckingham vista nell'introduzione. Sem- 
pre secondo la teoria, si ha che ci sono allora 5 — 3 = 2 paramentri adimen- 
sionali 



ni = -zh, n 2 = — , (14) 



EJS _ A 

pv z A A 

e quindi, per il Teorema II, la relazione funzionale generale che lega E c j S a 
A, p, v e A e 

dove $ e una funzione adimensionale indeterminata, per ottenere la quale e 
necessaria una vera teoria o un confronto diretto con i dati sperimentali. 

Ma qual'e il comportamento della funzione A<&(\/A) al variare di A e 
A? Euristicamente ci aspettiamo che se A e/o A aumentano, deve aumentare 
anche %. Infatti, se A e/o A aumentano a parita di v e p, le dimensioni fisiche 
dell'onda aumentano e con loro la sua massa e dunque la sua energia cinetica. 
Quindi, la funzione $(A/A) deve essere proporzionale al rapporto X/A per 
rendere conto dell'aumento di $ all'aumentare di A (con A fissato). Inoltre, 
AQ(\/A) deve crescere all'aumentare di A, con A costante. Ne consegue che 
all'aumentare di A, $(A/A) deve diminuire (e diminuisce poiche il rapporto 
X/A diminuisce) piu lentamente di A per ogni valore di A. Un semplice 
esempio di funzione $ che soddisfa queste caratteristiche potrebbe essere il 
seguente 

^ = k c pv 2 A a \^ (16) 
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dove k c e una costante adimensionale sconosciuta, a e (3 sono due numeri 
qualunque ma maggiori di zero e tali che a + f3 — 1. 

Come anticipato poco fa, l'onda marina trasporta anche energia poten- 
ziale gravitazionale, per il semplice fatto che una certa massa d'acqua, il picco 
dell'onda, e sollevata rispetto al pelo libero della superficie del mare. Cer- 
chiamo allora di derivare euristicamente l'energia potenziale dell'onda che 
attraversa l'unita di superficie, E p /S. Come prima, le quantita fisiche da 
cui Ep/S puo dipendere sono la densita p dell'acqua, che mutua la massa 
dell'onda, l'ampiezza A e la lunghezza d'onda A che esprimono la grandezza 
(altezza e lunghezza) dell'onda con cui abbiamo a che fare e infine la costante 
di acceleraione gravitazionale g, che tiene conto della forza di gravita rispetto 
alia quale l'onda possiede un'energia potenziale. 

Passando alle unita di misura (dimensioni), si ha: 

[E p /S] = M ■ T~ 2 , 

[p] = M-L~\ 

[A] = L, (17) 

[A] = L, 

[g] = l-t-\ 

Anche qui, scegliendo come parametri indipendenti A, p e g, si ha che ci 
sono 5 — 3 = 2 paramentri adimensionali 

n ' = ^- n -i- < 18 » 

e quindi, sempre secondo il Teorema IT, la relazione funzionale generate che 
lega Ep/S a A, p, g e A e 

dove \1/ e una funzione adimensionale indeterminata, per ottenere la quale, 
di nuovo, e necessaria una vera teoria o un confronto diretto con i dati speri- 
mentali. Anche in questo caso, il comportamento della funzione A 2 ^/(\/A) e 
molto simile a quello della funzione AQ(\/A). Infatti, al crescere di A e/o A 
l'onda cresce in dimensioni e piu l'onda e grande, maggiore e la sua energia 
potenziale, cioe 

Quindi, come prima, la funzione ^/(\/A) deve essere proporzionale al 
rapporto \/A per rendere conto dell'aumento di ^ all'aumentare di A (con 
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A fissato). Ma A 2 m(\/A) deve crescere all'aumentare di A, con A costante: 
ne consegue che all'aumentare di A, ty(\/A) deve diminuire piu lentamente 
di A 2 per ogni valore di A. Come prima, un semplice esempio di funzione $ 
che soddisfa queste caratteristiche potrebbe essere il seguente 

^ = k.pgA^X 5 , (20) 

dove, ancora, k p e una costante adimensionale sconosciuta e 7 e 5 due numeri 
positivi tali che 7 + 5 = 2. 

Ma ora torniamo agli tsunami. Se trascuriamo le dissipazioni legate 
all'attrito interno dell'acqua e con la sabbia della riva poco prima dell'impatto 
dell'onda, possiamo considerare l'energia totale dell'onda che attraversa l'unita 
di superficie, E t /S, come una quantita quasi costante durante la fase di for- 
mazione dello tsunami, cioe 

Et E r E„ 

i = i + f^ K - < 21 > 

per cui la relazione funzionale che ci fornisce % e 

pv 2 At> (j\ + pgA 2 ^ (jj ~ K t . (22) 

A questo punto, nei pressi della spiaggia v deve necessariamente diminuire 
(si avvicina ad una zona dove non c'e piu acqua!), cosi come deve ridursi la 
lunghezza d'onda A: in un certo senso l'onda marina rallenta e si comprime 
in prossimita della riva. Dall'equazione (j22j) risulta allora chiaro che se v 
e A diminuiscono, affinche E t /S rimanga quasi costante e necessario che A 
aumenti! Non c'e altra possibility di soddisfare tale vincolo, visti i com- 
portamenti di AQ(\/A) e A 2 ty(\/A) analizzati piu sopra; per ogni valore 
istantaneo di A, l'aumento di A implica l'aumento di A<f>(\/A) e A 2 ^(\/A). 

Per completezza va detto che la relazione (|22|) e stata ottenuta non 
tenendo conto dell'influenza che potrebbe avere la profondita del fondale 
(nelle equazioni (fTK|) e la profondita h non compare): questo ha reso 
la derivazione della (|2"2*|) alquanto diretta, ma ha anche confmato il range di 
applicability in quella che puo essere chiamata zona asintotica, quella zona 
dove il fondale non influisce poiche ha una profondita il cui valore numerico 
e molto piu grande di quello dell'ampiezza dell'onda (h ^> A). 

Ma allora dobbiamo credere o no alia dipendenza funzionale di A da v in 
prossimita della riva mostrata dalla (J22j) ? La risposta rigorosa e no, tuttavia 
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r Epicentro 




Figure 2: Esemplificazione grafica del fattore di attenuazione geometrica 
dell'energia dell'onda (eq. E3J). 

se A e originariamente (ovvero quando h ^> A) piu grande della lunghezza 
del tratto di mare in cui il fondale passa da h A a h ~ A, allora possiamo 
dire che la relazione ()22|) e accettabile. 

In ogni caso, la ()22|) ci mostra indicativamente qual'e la tendenza di A in 
funzione di v quando v diminuisce in prossimita della riva. 

A questo punto puo essere utile fornire qualche informazione relativa a 
tsunami verificatisi realmente in passato [2] . Molti degli tsunami di cui si ha 
documentazione storica (anche fotografica, vedi fig. EJ) hanno avuto origine 
nell'Oceano Pacifico. Nell'Oceano Pacifico, dove la profondita media delle 
acque e circa di 4000 metri, uno tsunami viaggia ad una velocita compresa 
tra i 700 e i 1000 chilometri orari (la velocita di crociera di un jet di linea). 
Inoltre, la sua ampiezza e spesso minore di un metro, mentre la sua lunghezza 
d'onda e mediamente di 100 chilometri (il passaggio di un'onda di tsunami 
in alto mare non puo essere praticamente percepita). In prossimita della 
riva, tuttavia, l'ampiezza dell'onda sale mediamente fino ad un massimo di 
30 metri! 

Passiamo ora a occuparci brevemente della quantita K t , l'energia to- 
tale dell'onda marina per unita di superficie attraversata. Supponiamo che 
un'onda anomala sia generata in alto mare da un evento naturale di grandi 
proporzioni (terremoto, impatto asteroidale, improvvise modificazioni oro- 
grafiche o eruzioni vulcaniche subacquee, etc.) e che l'energia trasferita al 
moto dell'acqua sia E (che sara solo una frazione dell'energia totale sprigion- 
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ata dall'evento). Ora, se il tratto di mare in cui si verifica il passaggio da 
li>Aa/i~A(ein cui c'e la riva) dista r dall'epicentro del fenomeno, 
allora possiamo supporre che K t abbia il seguente andamento funzionale 

K < K i (23) 

dove il fattore rA ^ tiene conto dell'attenuazione dovuta alia geometria: 
l'energia si distribuisce approssimativamente sul bordo di una circonferenza 
di raggio r e per un'altezza dell'ordine di A^>a (che e l'ampiezza dell'onda 
anomala prima che incontri la riva e si trasformi in tsunami), e quindi la 
superficie su cui si distribuisce e in qualche modo proporzionale a rA h ^A- 
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